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I. INTRODUCTION 

In our previous paper jl| (cited below as I) we carried out an accurate detailed study of 
the ground state l g of the molecular ion H 2 placed in a constant uniform magnetic field 
ranging from zero up to 4.414 x 10 13 G for all inclinations 0° — 90°. The goal of that study 
was to investigate the domain of existence of the H 2 ion. We showed that for all magnetic 
fields studied the molecular ion H 2 exists for moderate (not very large) deviations of the 
molecular axis from the magnetic field direction (moderate inclinations). Furthermore it was 
found that for each magnetic field the most stable configuration of minimum total energy 
corresponded to zero inclination, where the molecular axis coincides with magnetic field 
direction. We called this configuration the 'parallel configuration'. To this configuration the 
standard spectroscopic notation la g can be assigned. A major feature of this configuration is 
that with magnetic field growth the system becomes more and more bound (binding energy 
grows) and more and more compact (equilibrium distance decreases). 

The aim of the present paper is to continue the study initiated in I and to explore several 
low-lying excited states. At first we re-examine the ground state for the parallel configuration 
la g in the region 10 9 — 4.414 x 10 13 G. A detailed study of the la u state which is anti-bonding 
without a magnetic field is presented. Then the lowest states of different magnetic quantum 
numbers are investigated as well as the 2a g state. 

Atomic units are used throughout (h=m e =e=l) , although energies are expressed in Ry- 
dbergs (Ry). The magnetic field B is given in a.u. with Bq = 2.35 x 10 9 G [18]. 



II. GENERALITIES 



The Hamiltonian which describes two infinitely heavy protons and one electron placed in 
a uniform constant magnetic field directed along the z— axis, B = (0, 0, B) is given by (see 
e-g- Q) 

7i = p 2 + |- - — - — + ipA + Ap) + A 2 , (1) 

R r\ r 2 

(see Fig.l for notations), where p = —iV is the momentum, A is a vector potential which 
corresponds to the magnetic field B and is chosen in the symmetric gauge to be 

A=j(-y, x, 0) . 
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Hence the total energy Et of is defined as the total electronic energy plus the Coulomb 
energy of proton repulsion. In turn, the binding energy is defined as an affinity to having the 
electron as well as both protons being infinitely separated, E^ = B — Et- The dissociation 
energy is defined as an affinity to having one proton at infinity, E& = En — Et, where Eh 
is the total energy of the hydrogen atom in a magnetic field B. Spin degrees of freedom can 
be separated out and their analysis is straightforward. 

The problem is characterized by two integrals of motion: (i) angular momentum projec- 
tion on the magnetic field direction (z-direction) and (ii) spatial parity p. Sometimes the 
parity a corresponding to interchange of charged centers 1 <-> 2 is used, which is connected 
with the magnetic quantum number and spatial parity, 

p = o-(-L)H . 

If the case m is even, both parities coincide, p = a. Thus, any eigenstate has two definite 
quantum numbers: the magnetic quantum number m and the parity p with respect z — > 
—z. Therefore the space of eigenstates is split into subspaces (sectors) each of them is 
characterized by definite m and a. Notation for the state we are going to use is based on the 
following convention: the first number corresponds to the number of excitation - "principal 
quantum number", e.g. the number 1 is assigned to the ground state, then a Greek letter 
a, 7i,5 corresponds to m — 0,-1,-2, respectively, with subscript g/u (gerade/ungerade) 
describing positive/negative parity with respect z — > —z. 

Most of the excited states we study are the lowest states (of the type of the ground state) 
of the sectors with different magnetic quantum numbers m and p. It is quite obvious from 
the physical point of view that the ground states of the sectors with m > always have 
larger total energies than those with m < 0. Therefore we restrict our consideration to the 
states with m — 0, — 1, —2. 

Conservation of the z-component of the angular momentum assumes the wave function 
of the electron (in cylindrical coordinates (p, <p, z)) can be taken in the representation: 

where m is magnetic quantum number. Let us gauge rotate the Hamiltonian 

7i m = e - im V~ H fte im V H =P 2 m + ^----+mB + - , (2) 

R r\ r 2 4 
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-R/2 R/2 * 
► 

FIG. 1: Geometrical setting for the ion placed in a magnetic field directed along the z-axis. 
The protons are situated in the y — z plane at a distance R from each other and marked by bullets. 



where 



p m = e -vtn>p p -M p e *m> p \m\ 



is the gauge rotated momentum (covariant momentum). The constant term mB describes 
the linear Zeeman effect splitting. It can be absorbed to a definition of total energy. The 
representation (J2J) is rather convenient since each Hamiltonian for fixed m describes the 
family of eigenstates with quantum number m and can be treated independently of the states 
with wl different from m. Now the Hamiltonian (J2J) has only the invariance corresponding 
to the spatial parity conservation. 

We are going to use the variational method in a way similar to what was done in I. The 
recipe of choice of trial function is based on physical arguments and is described in full 
generality in P (see for details the article I). Eventually, the ground state trial function for 
fixed m and a is chosen in a form 

^(trial) = A ^ + + ? (3) 



4 



where 



e -ai(n+nO e -B/?ip» ^ ifa = +l, 
0, if er = — 1 , 

and cr = ±1, m = 0, ±1, ±2 . . .. Here ^1,2,3 and 01,2,3,4, @\. 2,3 as well as i? are variational 
parameters, which are certainly different for different m [19(. The functions ^1^,3 carry a 
certain physical meaning. They describe coherent (incoherent) interaction of the electron 
with the protons as well as their non-linear interpolation, respectively. Calculations were 
performed using the minimization package MINUIT from CERN-LIB. Numerical integra- 
tions were carried out with a relative accuracy of ~ 10 -9 by use of the adaptive NAG-LIB 
(D01FCF) routine. All calculations were performed on dual PC Pentium-4, with two pro- 
cessors of 2.8 GHz each. Every particular calculation of given eigenstate at fixed magnetic 
field including minimization has taken in total about an hour of CPU time. However, when 
the variational parameters are found it take a few seconds of CPU time to calculate the 
variational energy. 

It is necessary to mention two technical difficulties we encountered. Calculation of two- 
dimensional integrals with high accuracy which appeared in the problem has required a 
development of a very sophisticated numerical technique. We created a 'dynamical par- 
titioning' of the domain of integration, which depend on values of variational parameters 
similar to what was done in I. The domain partitioning was changed with a change of the pa- 
rameters. Sometimes the number of sub-domains was around 50. Another technical problem 
is related with very complicated profile of variational energy as the function of variational 
parameters which is characterized by many local minima, saddle points and valleys. Local- 
ization of the global minimum numerically of such a complicated function with high accuracy 
is difficult technical problem which becomes even more difficult in the case of ten or more 
variational parameters. Examining the physical relevance of trial functions allows one to 
avoid spurious minima. The parameters obtained in (J3J) at every step of minimization were 
always examined from the physical point of view. Such considerations are always something 
of an art. 
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III. m = 



The m = subspace consists of two subspaces, a = 1 (even states) and a = — 1 (odd 
states). 

A. la g state (a = 1) 

The state l g was thoroughly investigated in the paper I for the whole range of inclinations 
= 0° — 90° (for settings see below Fig.2). At 9 = 0° this state becomes the state la g and 
our variational anzatz ^\ ag © describing this state depends on ten parameters. As was 
mentioned above, the search for the global minimum numerically with high accuracy in the 
case of so many variational parameters is a difficult technical task. Although this state was 
thoroughly studied in [4] we decided to repeat the calculations using a more sophisticated 
strategy for localizing the minimum. The essential new element of the strategy was to im- 
pose an extra (natural) condition that the variational parameters change smoothly with B. 
Finally, it led to an improvement of the results in comparison to I and to previous calcu- 
lations. It is worth mentioning that this recalculation is very important for calculation of 
the excited 2a g state, where the orthogonality condition on trial functions must be imposed, 
(ipia , ip2a g ) = 0. It is evident that an intrinsic inaccuracy in ip 1(Tg due to the approximate 
nature of the trial function JHJ) as a function of x is a source of inaccuracy in the energy 
of the ?/>2o- 9 state. Thus, a reduction of this inaccuracy requires knowledge of the function 
ipi ag as accurately as possible. The above-mentioned strategy allowed us to improve our 
previous results reported in [4] on total and binding energies (see Table I) and also on lowest 



rotational- vibrational energies (see Table II). Qualitative conclusions obtained in 
unchanged. 



4] 



remain 
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TABLE I: Total Ex and binding energies, and equilibrium distance R eq for the state 1<7 9 , which 
is the global ground state of the H2 ion. Error bars for the equilibrium distance indicate a domain 
in R where the value of binding energy remains the same within the indicated number of their 
digits shown in the table. 
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TABLE II: Energies of the lowest vibrational (E v n,) and rotational (E rot ) electronic states associ- 
ated with the lag state. 
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B. 2a g state (<r = 1) 



This is the first excited state in the family of states with quantum numbers m = 0, a = 1. 
Here the trial function is taken in the form 



^2 



A^x + A 2 ^ 2 + A 3 i/j 3 



(4) 



with 



ipi = (n + r 2 -C x )e 



-ai(n+r 2 )-/3iSp 2 /4 



1p2 



[n - C 2 )e~ & ^ + (r 2 - C 2 )e 



OL2T-2 



-03 



-/3 2 Bp 2 /4 



-OLZT1— 04TI 



-PsBp 2 /4 



(n + ar 2 - C 3 )e- a3ri - Q4r2 + (r 2 + an - C 3 )e" 

(cf. ©), where v4 lj2 3 and ai^A, $1,2,3, 0, £-1,2,3 as well as i? are variational parameters 
This eigenfunction should be orthogonal to the ipi a state trial function found in the previous 
Section. The total number of variational parameters in (jlj) is 13. 

The results obtained are presented in Table III. This state is characterized by much 
smaller binding energy compared to the la g state and is much more extended. The binding 
energy displays a rather slow increase while the equilibrium distance decreases slowly as the 
magnetic field grows. This excited state is unstable with respect to dissociation to H + p. 
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TABLE III: Total and binding energies, and equilibrium distance R eq for the state 2a g 
( m = 0,a = 1). 
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C. la u state (<r = — 1) 



In the absence of a magnetic field the la u state (m = 0, a — —1) is essentially repulsive 
and antibonding. However, in a strong magnetic field this state becomes bound. Due to 
this fact we want to study this state in full generality, for different magnetic fields and 
inclinations. 

In the absence of a magnetic field, the \o u state is characterized by a shallow minimum 
in the total energy situated at large internuclear distance (see, for example, [3] (a), 
Also this state is a weakly bound state with respect to dissociation and it becomes even 
unbound if nuclear motion is taken into account. So far not many studies have been carried 
out for this state. Our major finding is that in the presence of a magnetic field of the 
magnitude 10 9 < B < 4.414 x 10 13 G the total energy surface of the system (ppe) in the 
state la u exhibits a well-pronounced minimum. Similar to the la g state, both total (Et) 
and binding (Eb) energies of the la u state increase as the magnetic field grows, while the 
equilibrium distance decreases. However, the accuracy of our calculations does not allow us 
to make a definitive conclusion about the stability of the system in this state with respect 
to dissociation and nuclear motion effects. In the case of non-zero inclination 6^0° (for 
definition see Fig. 2) we denote this state as l u reflecting the fact that the only parity 
conservation exists. In I it was shown that for B > 10 11 G and large inclinations the l g state 
disappears and hence the molecular ion H% does not exist. Thus, it seems it makes no sense 
to study the l u state in this domain. We checked a consistency of this statement verifying 
that always inequality E^ 9 (R) < E^ CTu (R) holds. 

To study the l u state we use the following form of the vector potential corresponding to 
a constant magnetic field B = (0, 0, B) 

A=B((£-l)y, £r, 0) , (5) 

where £ is a parameter, which later will be considered as variational. If £ = 1/2 we get the 
well-known and widely used gauge which is called symmetric or circular. If £ = or 1, we 
get the asymmetric or Landau gauge (see Q]). By substituting © into (JTJ) we arrive at a 
Hamiltonian of the form 



H = -V 2 + | - - - - - 2iB[(£ - l)yd x + &d y ] + B 2 [ex 2 + (1 - 0¥] • (6) 
R r± r-i 
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FIG. 2: Geometrical setting for the inclined ion in a magnetic field directed along the z-axis. 



The trial function is chosen in the form 



u 



Ai^a + A 2 ^ 2 , 



(7) 



with 



ih. = (e 




e -B[/3 lx tix 2 +f3 ly (l-Oy 2 ] 



^2 = (e 



— «2f 1— 03^2 



— e 



— «2»*2 — 03^1 



B[/32^x 2 +/3 2 „(i-Oy 2 ] 



where Ai, A 2 are parameters and one of them is kept fixed by a normalization condition. All 
parameters 0:1,2,3, Pix,\y,2x,2yi A\, A 2 and £ are variational parameters. It is evident that if 
9 = 0°, the rotational invariance along z-axis exists and the vector potential should be taken 
in a form supporting this invariance. Hence the parameter £ in (jSJ) takes value £ = 1/2 and 
the parameters f3\ x = Piy,p 2 x = fiiy 

Numerical study for the l u state was carried out for different inclinations with the results 
at 0°, 45° and 90° for magnetic fields B = - 4.414 x 10 13 G as shown in Tables IV- VI. The 
immediate conclusion is that 



for all magnetic fields, where this comparison makes sense (see below). Hence, similar to the 
l g state, the highest molecular stability of the l u state occurs for the parallel configuration, 
at 9 = 0° (see I). Also, the binding energy growth is maximal as a function of magnetic field 
for the parallel configuration. Therefore, the stability of H 2 in the parallel configuration 
in the l u state increases as the magnetic field grows, again similarly to what happens for 
the l g state. These results suggest the following picture for appearance of a bound state for 



E T (0°) < £ T (45°) < E T (90°) 
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TABLE IV: la u state in the parallel configuration, 9 = 0°. Total (Et) and binding (Eb) energies 
are in Ry and equilibrium distance R eq in a.u. 



B 



E T (Ry) 



E b (Ry) 



Reg (a.U.) 



B = 

10 9 G 

1 a.u. 

10 10 G 

10 a.u. 
10 n G 
100 a.u. 
10 12 G 
1000 a.u. 
10 13 G 
10000 a.u. 
4.414 x 10 13 G 



-1.00010 
-1.00012 

-0.92103 
-0.917134 

-0.66271 
-0.66 

1.63989 
2.1294 

6.52362 

36.8367 

92.4257 

413.6175 

984.6852 

4232.554 

9971.727 

18750.07 



1.00010 
1.00012 

1.34656 

1.66271 

I. 66 

2.61500 

3.47638 

5.7165 

7.5743 

II. 9144 
15.3148 
22.765 
28.273 
32.912 



12.746 
12.55 

11.19 
10.55 

9.73 
9.6 

7.18 
4.18 

6.336 

4.629 

3.976 

3.209 

2.862 

2.360 

2.134 

2.021 



Lopez et al [4] 
Peek-Katriel [14] 

Present 

Peek-Katriel [14] 
Present 

Kappes et al [7] 
Present 

Peek-Katriel [14] 

Present 

Present 

Present 

Present 

Present 

Present 

Present 

Present 



the l u state: for small magnetic fields the minimum in the total energy arises at very large 
internuclear distances j^, then, as the magnetic field grows, the position of the minimum 
moves to smaller and smaller internuclear distances. 

Our results for B > and 8 = 0° give the lowest total energies compared to other 
calculations. In general, they are in a good agreement with those by Kappes-Schmelcher 
Q (a) as well as by Peek-Katriel Q for B = 0, 10 9 G, although for B = 10 10 G a certain 
disagreement is observed (see Table IV). However, for 9 = 90° our results are in striking, 
qualitative contrast with those by Wille 5], where even the optimal configuration is attached 
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TABLE V: ' 


rotal Ex, binding 


Eb energ 


ies and equilibrium distance R eq for the l u i 


state in the 


configuration 


9 = 45°. Optimal value for the gauge parameter £ 


is shown (see 


text). 




B 


E T (Ry) 




E b (Ry) 


Req (a.U.) 






10 9 G 


-0 870391 




1 295923 


8.053 




0.9308 


1 a.u. 


-0 509041 




1 509041 


6.587 




0.9406 


10 10 G 


2 267998 




1 987321 


4.812 




0.9671 


10 a.u. 


7.692812 




2.307188 


4.196 




0.9808 


10 11 G 


39.71061 




2.84258 


3 538 




0.9935 


100 a.u. 


96.88464 




3.11536 


3.278 




0.9968 


10 12 G 


422.0074 




3.5245 


3.020 




0.9991 


1000 a.u. 


996.3044 




3.6956 


2.894 




0.9996 


TABLE VI: Total Et, binding Ef, energies and equilibrium distance R eq for the l u state at 8 = 90°. 


Optimal value for the gauge parameter £ 


is shown (see text). 
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10 9 G 


-0.867234 


1.292766 


8.784 


0.9692 


Present 




1 a.u. 


-0.49963 


1.49963 


7.264 


0.9737 


Present 






-0.65998 


1.65998 


5.45 




Kappes 


et al [^(b) 


10 10 G 


2.29365 


1.96167 


5.517 


0.9866 


Present 




10 a.u. 


7.72998 


2.27002 


4.872 


0.9923 


Present 




10 11 G 


39.76500 


2.78819 


4.154 


0.9975 


Present 




100 a.u. 


96.93497 


3.06503 


3.875 


0.9988 


Present 





to 9 = 90°, contrary to our conclusion. For instance, at B — 10 10 G in j| the values 
Et, = 2.593 Ry and R eq = 2.284 a.u. are given, while our results are = 1.9617 Ry 
and R eq = 5.517 a.u., respectively (see Table VI). Similar, but less drastic disagreement, 
is observed with the results in p)®. We can only guess this disagreement is due to the 
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shallow nature of the minimum, but a real explanation of this fact is missing. Independent 
calculations are needed in order to resolve this contradiction. 

The analysis of Tables IV- VI shows that for 9 > 0° and fixed magnetic field the total 
energy of in the \ u state is always larger than the total energy of the hydrogen atom 
. It means that the if^-ion in the l u state is unstable towards dissociation to H + p. 
For 9 ~ 0° the total energies presented for the H2 ion and the most accurate results for the 



hydrogen atom 16[ are comparable within the order of magnitude 10~ 4 — 10 -5 . We estimate 
that the accuracy of our calculations is of the same order of magnitude — 10 -5 . This 
prevents us from making a conclusion about the stability of in the l u state with respect 
to dissociation. Thus, the only reliable conclusion can be drawn that the minimum is very 
shallow. 

The l u state is much more extended than the l g state: for fixed magnetic field the 
equilibrium distance of the l g state is much smaller than that for the l u state. This picture 
remains the same for any inclination. It is quite impressive to observe the much lower rate 
of decrease of R eq in the range B = — 4.414 x 10 13 G with magnetic field growth. For 
example, in the case of the parallel configuration, 9 = 0°, for the state l u the equilibrium 
distance falls ~ 6 times compared to the l g state, where it falls ~ 20 times. 

The behavior of the equilibrium distance R eq of the l u state as a function of inclination 
is quite non-trivial (see Tables IV- VI). As in the l g state, the H^~-ion in the l u state for 
B < 10 12 G is most extended in the parallel configuration. 

IV. m = -1 

The subspace consists of two subspaces, o = 1 (even states) and o = — 1 (odd states). 
A. lir u state (<r = 1) 

In order to study the l7r M (m = —1 and a — 1) state we take the trial function Q. 
The results are presented in Table VII. In general, our results are more accurate than those 
obtained in other calculations giving lower total (and correspondingly, the higher binding) 
energies with the onlv exception of the magnetic field 5 = 10 a.u. where the result for 
binding energy from [5J] is better in the fourth digit. The results for B = 1 a.u. obtained 
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TABLE VII: Total Et, binding energies and equilibrium distance R eq for the excited state 1-7T, 
(m = —1, a = 1). 



B 


Et (By) 


E b (Ry) 


R eq (&.u.) 




10 9 G 


-0.293592 


0.719123 


4.940 


Present 


1 a.u. 


-0.020150 


1.020150 


3.676 


Present 




-0.02014 


- 


3.68 


Kappes et al [7] 




-0.02011 


- 


3.75 


Wille [5] 


10 10 G 


2.371845 


1.883474 


2.130 


Present 


10 a.u. 


7.29682 


2.70318 


1.526 


Present 






2.6862 


1.510 


Vincke-Baye [8] 






2.7046 


1.510 


Wille [5J 


10 11 G 


37.6490 


4.9042 


0.887 


Present 


100 a.u. 


93.1127 


6.8873 


0.651 


Present 






6.8774 


0.645 


Vincke-Baye [8] 






6.8548 


0.645 


Wille [5] 


10 12 G 


413.6306 


11.902 


0.395 


Present 


1000 a.u. 


983.874 


16.126 


0.301 


Present 


10 13 G 


4229.183 


26.136 


0.195 


Present 


10000 a.u. 


9965.932 


34.068 


0.154 


Present 


4.414 x 10 13 G 


18741.89 


41.09 


0.130 


Present 



in pj are not very precise in R eq (see Table I therein), which explain their difference with 
the results by others. The binding energy at B — 10 and 100 a.u. in is calculated 
for the same equilibrium distances as those found in Q| (see Table IV in Ref. j^). Like 
for all studied states the binding energy grows steadily with magnetic field increase while 
the equilibrium distance shrinks in a quite drastic manner. If for small magnetic fields the 
equilibrium distance R eq is several times larger than this distance for the la g state, for large 
magnetic fields these equilibrium distances become comparable. Among m = — 1 states the 
state 1tt u has the smallest total energy. 
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TABLE VIII: Total Et, binding Ef, energies and equilibrium distance R eq for the state 1-K g (m = 
-l,a = -l). 



D 


E T (Ry) 


E b (Ry) 


R eq (a.u.) 




10 9 G 


-0.232060 


0.65759 


20.10 


Present 


1 a.u. 


0.086868 


0.91313 


14.05 


Present 




0.0866 




13.5 


Kappes et al [7J 


10 10 G 


2.641122 


1.61420 


9.370 


Present 


10 a.u. 


7.749819 


2.25018 


7.622 


Present 


10 11 G 


38.67642 


3.87677 


5.622 


Present 


100 a.u. 


94.73386 


5.26614 


4.791 


Present 


10 12 G 


416.9354 


8.59654 


3.767 


Present 


1000 a.u. 


988.7286 


11.2714 


3.321 


Present 


10 13 G 


4238.038 


17.2810 


2.708 


Present 


10000 a.u. 


9978.175 


21.8254 


2.420 


Present 


4.414 x 10 13 G 


18757.273 


25.7054 


2.237 


Present 



B. 



lng state (<t = — 1) 



In order to study the ln g state (m = — 1 and a = —1) we take the trial function 
results are presented in Table VIII. For B = 1 a.u. our total energy deviates from 
third digit and an independent calculation would be desirable. 



j). The 

in the 
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V. m = -2 



The subspace consists of two subspaces, a = 1 (even states) and a = — 1 (odd states). 
A. 15 g state (a = 1) 

In order to study the 15 g state (m = —2 and a = 1) we take the trial function @. The 
results are presented in Table IX. In |5] for B = 1 a.u. the equilibrium distance is simply 
placed equal to 5.0 a.u. (see Table I therein). For B = 10, 100 a.u. the energies computed 
in (i| were calculated for the same equilibrium distances as those found in [^] (see Table 
IV in (5I). Among m = —2 states the 15 g state has the smallest total energy. It is worth 
mentioning a drastic decrease of R eq with magnetic field growth similar to what appears for 
la g and 1tt u states. 
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IAdHj IX: lotai hiTt 


binding E b enerj 


$ies and equilibrium 


distance R eq 


for the state 15g (m = 


—2, a — +L). 










D 


E T (Ry) 


E b (Ry) 


.R e9 (a.u.) 




10 9 G 


-0.107945 


0.533477 


6.865 


Present 


1 a.u. 


0.221163 


0.778837 


4.872 


Present 




0.22112 




4.87 


Kappes et al [7] 




0.22126 




5.0 


Wille [5J 


10 10 G 


2.77538 


1.47994 


2.694 


Present 


10 a.u. 


7.85113 


2.14887 


1.907 


Present 






2.1306 


1.880 


Vincke-Baye [8] 






2.1496 


1.880 


Wille [5J 


10 11 G 


38.58470 


3.9685 


1.080 


Present 


100 a.u. 


94.38093 


5.6191 


0.782 


Present 






5.6058 


0.778 


Vincke-Baye [8] 






5.510 


0.778 


Wille [5] 


10 12 G 


415.6710 


9.8609 


0.470 


Present 


1UUU a.u. 


986.5119 


13.4881 


U.353 


Present 


10 13 G 


4233.125 


22.194 


0.225 


Present 


10000 a.u. 


9970.802 


29.198 


0.176 


Present 


4.414 x 10 13 G 


18747.572 


35.407 


0.148 


Present 
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TABLE X: Total Et, binding E\, energies and equilibrium distance R eq for the state 15 u (m = 



-2,a = -l). 



JJ 


Et (fty) 




-Keg (a.U. J 




1U (or 


n n^Q7Q 
-U.Uoo i o 


n /i n/i of; 


oq ono 


Present 


1 a.u. 


n on /1 1 n 
0.zy410 


u. /ooyo 


lb. 3 / / 


Present 




n oni^ 
O.zyoo 




i n 
lo.O 


Kappes et al [7] 


11) b 


o n^T/io 

z.yv /42 


i.z r /yu 


11/1 

11.4/5 


Present 


11) a.u. 


o.iyoyz 


i cm nc 
1.80108 


n /i K q 
y.458 


Present 


i nil 

ID Lr 




3.14/23 


f; qc;q 


Present 


inn a n 


Q5 6Q549 


4 30458 


5 61 Q 




10 12 G 


418.4335 


7.0984 


4.071 


Present 


1000 a.u. 


990.6416 


9.3584 


3.406 


Present 


10 13 G 


4240.834 


14.485 


2.625 


Present 


10000 a.u. 


9981.587 


18.413 


2.391 


Present 


4.414 x 10 13 G 


18761.18 


21.80 


2.230 


Present 



B. 1S M state (er = — 1) 

In order to study the 15 U state (m = —2 and a = — 1) we take the trial function The 
results are presented in Table X. 
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VI. DISCUSSION 



In Table XI a summary of total energies of eigenstates explored in this article for magnetic 
fields ranging from 10 9 G to 4.414 x 10 13 G is presented. An analysis of Table IX allows to 
draw a certain immediate conclusions: 

1. The state la g is the global ground state for all magnetic fields. It is rather evident 
that this statement remains valid in general, when even the states other than studied 
are taken into account (Perron theorem); 

2. For the states with fixed m the lowest total energy corresponds to the state of positive 
parity a = +1. We guess that this statement remains correct in general; 

3. For the same parity a ground states are ordered following the value of m, 

E? 9 < < E l T &9 , 
E%? u < Et 9 < E T Su . 

4. There exist several true level crossings, 

• For B < 10 12 G 

Efi u < Eft" , 

while for B > 10 12 G 

E l T au > E^ u , 

• For B < 10000 a.u. 
while for B > 10000 a.u 

• For B < 10 a.u. 
while for B > 10 a.u 



E£ u < E? 
E x ^ > E% 



Ulrp < Ulrp , 



P l7r 9 \ J? 1S 3 

Ulrp > Ulrp , 
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• For B < 1 a.u. 

Hi j, ^_ Hi j, ? 





while for B > 1 


a.u. 


E l T &9 

1 


< E 2 ° a , 

1 ' 








• 


For B < 1 a.u. 






> E 2 t 9 , 










while for B > 1 a.u. 


E]? u 


< E 2 T ° a . 








TABLE XI: Comparison of the total energies Et (in 


Rydbergs) for the low-lying states of the 


molecular 


ion for magnetic fields 10 9 G - 


■ 4.414 x 


10 1 










B 


la g 


1ct u 


lvr„ 




lTTg 


is g 






10 9 G 


-1.15070 


-0.92103 


-0.29359 




-0.232060 


-0.107945 


-0.068727 


-0.121343 


1 a.u. 


-0.94992 


-0.66271 


-0.20150 




0.086868 


0.22117 


0.29410 


0.34912 


10 10 G 


1.09044 


1.63989 


2.371845 




2.641122 


2.77538 


2.977418 


3.39938 


10 a.u. 


5.65024 


6.52362 


7.296816 




7.749819 


7.85113 


8.198922 


9.02452 


10 11 G 


35.04320 


36.83671 


37.64895 




38.67642 


38.58470 


39.40596 


41.4090 


100 a.u. 


89.7090 


92.4257 


93.11267 




94.7339 


94.38093 


95.69542 


98.7822 


10 12 G 


408.3894 


413.6175 


413.6306 




416.9354 


415.6710 


418.4335 


424.2278 


1000 a.u. 


977.2214 


984.685 


983.874 




988.7286 


986.5119 


990.6416 


998.662 


10 13 G 


4219.565 


4232.554 


4229.183 




4238.038 


4233.126 


4240.834 


4253.937 


10000 a.u. 


9954.203 


9971.727 


9965.932 




9978.175 


9970.802 


9981.587 


9998.608 


4.414 X 10 13 


G 18728.477 


18750.070 


18741.889 




18757.273 


18747.572 


18761.180 


18781.576 
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VII. CONCLUSION 



We have carried out an accurate, non-relativistic calculation in the Born-Oppenheimer 
approximation for the low-lying states of the if^ 1- molecular ion in the parallel configuration 
at equilibrium in the framework of a unique computational approach. The lo~ u state is 
considered in full generality for all inclinations of the molecular axis vs. magnetic field 
direction. We studied constant uniform magnetic fields ranging from B = 10 9 G up to 
B = 4.414 x 10 13 G, where non-relativistic considerations hold, although our method can 
be naturally applied to study the domain B < 10 9 G. We used a variational method with 
a very simple trial function with a few variational parameters inspired by the underlying 
physics of the problem. Thus our trial function can be easily analyzed and in contrast to 
other approaches our results can be easily reproduced. The trial function (3) can be easily 
modified to explore other excited states. 

The present study of several low-lying excited states complements a study of the ground 
state performed in I. Usually the total, binding, dissociation and transition energies grow 
with increase in the magnetic field, reaching values of several hundred eV at magnetic fields 
of 10 12 — 10 13 G. These results can be used to construct a model of the atmosphere of an 
isolated neutron star 1E1207. 4-5209 (see This will be done elsewhere. 
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